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The order parameter for the pure Yang-Mills phase transition is the Polyakov loop which encodes 
the symmetries of the Zn center of the SU{N) gauge group. On the other side the physical degrees 
of freedom of any asymptotically free gauge theory are hadronic states. Using the Yang-Mills trace 
anomaly and the exact Zn symmetry we construct a model able to communicate to the hadrons 
the information carried by the order parameter. 
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I. INTRODUCTION 

Investigating the SU{N) deconfinement phase transi- 
tion is, in general, a complex problem. At zero quark 
density importance sampling lattice simulations are able 
to provide vital information about the nature of the tem- 
perature driven phase transition for 2 and 3 colors Yang- 
Mills theories with and without matter fields (see |M 2] 
for 3 colors) . Different a p proaches [M H, B[_H, 0, ^, £ , 



are used 
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in literature to tackle/study the features of this phase 
transition. 

At zero temperature SU{N) Yang-Mills theory is 
asymptotically free and the physical spectrum of the the- 
ory consists of a tower of hadronic states referred as glue- 
ball and pseudo-scalar glueballs. The theory develops a 
mass gap and the lightest glueball has a mass of the or- 
der of few times the confining scale. The classical theory 
is conformal while quantum corrections lead to a non- 
vanishing trace of the energy momentum tensor. 

At nonzero temperature the Zn center of SU{N) is a 
relevant global symmetry |^ and it is possible to con- 
struct a number of gauge invariant operators charged 
under Zn among which the most notable one is the 
Polyakov loop: 
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V denotes path ordering, g is the SU{N) coupling con- 
stant and X is the coordinate for the three spatial di- 
mensions while T is the euclidean time. The i field is 
real for N — 2 while otherwise complex. This object 
is charged with respect to the center Zn of the SU (N) 
gauge group |22 under which it transforms as £ ^ zi 
with z 



e z^ 
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A relevant feature of the Polyakov loop 
is that its expectation value vanishes in the low temper- 
ature regime while is non zero in the high temperature 
phase. The Polyakov loop is a suitable order parameter 
for the Yang-Mills temperature driven phase transition 

H- 

This behavior recently led Pisarski to model the 
Yang-Mills phase transition as a mean field theory of 
Polyakov loops. This model is often referred as the 
Polyakov Loop Model (PLM). Using this model one can 
infer that the SU{2) phase transition is second order 



while a phase transition (as function of the tempera- 
ture) is a weak first order for SU{3). The predictions 
are in reasonable agreement with lattice results. More- 
over the PLM is used to model the Yang-Mills free en- 
ergy. Recently some interesting phenomenological PLM 
inspired models aimed to understand RHIC physics were 
constructed ^ 

Here we will consider pure gluon dynamics. This al- 
lows us to have a well defined framework where the Zn 
symmetry is exact. The hadronic states are the glue balls 
fields (H) and their effective theory at the tree level is 
constrained by the Yang-Mills trace anomaly. 

A real puzzle to me is how the information about 
the Yang-Mills phase transition encoded, for example, 
in the Zn global symmetry can be communicated to the 
hadronic states of the theory. Here we propose a concrete 
model which can help resolving this puzzle. 

This model is constructed using trace anomaly and the 
Zn symmetry. We will demonstrate that the informa- 
tion carried by £ is efficiently transferred to the glueballs. 
More generally the glueball field is a function of £: 



H = HI 



(2) 



Our results can be tested via first principle lattice sim- 
ulations p3| and support the recent phenomenological 
investigations_ll9, |H 
In section 



we present the model. In section |l^ we 
consider the two colors Yang-Mills theory while in IV the 
three color theory is considered. We finally conclude in 
section ^ 



II. THE MODEL 

The hadronic states of the Yang-Mills theory are the 
glueballs. At zero temperature the Yang-Mills trace 
anomaly has been used to constrain the potential of the 
lightest gluebaU state H p3| : 
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A is chosen to be the confining scale of the theory and 
H is a mass dimension four field. This potential cor- 
rectly saturates the trace anomaly when H is assumed 
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to be proportional to Tr [Gf^„G^'^] and G^^i, is the stan- 
dard Yang-Mills field strength. The potential nicely en- 
codes the properties of the Yang-Mills vacuum at zero 
temperature and it has been used to deduce a number of 
phenomenological results . 

At high temperature Pisarski conjectured that the 
Yang-Mills pressure can be written in terms of the field 
£. This free energy must be invariant under Zjv and it 
takes the general form: 



V[£] ^ T'^T[£]. 



(4) 



is a polynomial in £ invariant under Zn and the co- 
efficients depend on the temperature itself allowing for 
a mean field description of the Yang- Mills phase transi- 
tions. 

We now marry the two models by requiring both fields 
to be present simultaneously at non zero temperature. 
The theory must reproduce at zero and low tempera- 
tures the ordinary glueball Lagrangian while the PLM 
at high temperatures. We propose the following effective 
potential: 



V[H,£] = ^ln 



H 

A4 



Vt[H] + HV[£]+T^V[£], (5) 



where V [£] and V [£] are general (but real) polynomials 
in £ invariant under whose coefficients depend on the 
temperature. The explicit dependence is not known and 
should be fit to lattice data. Dimensional analysis and 
analyticity in H when coupling it with £ severely restricts 
the effective potential terms. We stress that HV[£] is the 
most general interaction term which can be constructed 
without spoiling the zero temperature trace anomaly. 

Further nonanalytic interaction terms can arise when 
considering thermal and quantum corrections and are 
partially contained in Vt[H] which schematically rep- 
resents the temperature of a gas of glueballs. In the 
following we will not investigate in detail such a term. 
Our theory cannot be the full story since we neglected 
(as customary) all of the tower of glueballs and pseudo- 
scalar glueballs as well as the infinite series of dimension- 
less gauge invariant operators with different charges with 
respect to Z^. Nevertheless the potential is sufficiently 
general to hope to capture the essential features of the 
Yang-Mills phase transition. 

When the temperature T is much less than the con- 
fining scale A the last term in Eq. (H) can be safely ne- 
glected. Since the glueballs are relatively heavy com- 
pared to the A scale their temperature contribution 
VT[iJ] can also be disregarded. At low temperatures the 
theory reduces to the standard glueball potential aug- 
mented by the third term which does not affect the trace 
anomaly. 

At very high temperature (compared to A) the last 
term dominates {H itself is very small) recovering the 
picture in which £ dominates the free energy. In this 
regime we have J-'[£] = V[£]. 



We can, in principle, compute all of the relevant ther- 
modynamical quantities in our approach, i.e. entropy, 
pressure etc. 

A relevant object is the trace of the energy-momentum 
tensor 0^. At zero temperature and when the poten- 
tial is a general function of a set of bosonic fields 
with mass-dimensions d„ one can construct the associ- 
ated trace of the energy-momentum tensor via: 



9^ 



At finite temperature we still define our temperature de- 
pendent energy- momentum tensor as in Eq. (|^). Here 
H possesses engineering mass dimensions 4 while £ is 
dimensionless yielding the following temperature depen- 
dent stress energy tensor: 



e'f.iT) = -2H + AT'^V[£] + 4 



1 -H- 
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m 



Vt[H]. (7) 



is normalized such that (0|8^|0) = e — 3p with e the 
vacuum energy density and p the pressure. At zero tem- 
perature only the first term survives yielding magnetic 
type condensation typical of a confining phase while at 
extremely high temperature the second term dominates 
displaying an energy density and pressure of the decon- 
fined phase. 

The theory containing just £ can be obtained integrat- 
ing out H via the equation of motion: 



5V[H, £] 
5H 



0. 



(8) 



Formally this is justifiable if the glueballs degrees of free- 
dom are very heavy. For simplicity we neglect the contri- 
bution of Vt [H] as well as the mean- field theory correc- 
tions for £. However in the future a more careful treat- 
ment which also includes the kinetic terms should be 
considered. Within these approximations the equation 
of motion yields: 



A" 

H[£] = — exp[-2P[ 
e 



(9) 



The previous expression shows the intimate relation be- 
tween £ and the physical states of strongly interacting 
theories. 

After substituting Eq. (0) back into the potential (|^) 
and having neglected Vt[^ we have: 



A"* 

V[£] = T'^V[£] - — exp [-2r 
2e 



(10) 



This formula shows that for large temperatures the only 
relevant energy scale is T and we recover the PLM model. 
However at low temperatures the scale A allows for new 
terms in the Lagrangian. Besides the and the A* 
terms we also expect terms with coefficients of the type 
TA^ and T^A^ and T'^A. However in our simple model 
these terms do not seem to emerge. 
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Expanding the exponential we have: 



V[l] = T^V[l] + —V[l] 



A4 
2e 



(11) 



Since V[P\ and V[P\ are real polynomials in £ invariant 
under Zj^ we immediately recover a general potential in 



III. THE TWO COLOR THEORY 

To illustrate how our formalism works we first consider 
in detail the case N = 2 and neglect for simplicity the 
term Vr [H] . This theory has been extensively studied via 
lattice simulations |2^, and it constitutes the natural 
playground to test our model. Here £ is a real field and 
the Z2 invariant V[£] and P[i] are taken to be: 



(12) 



with fli, 02 and bi unknown temperature dependent 
functions which should be derived directly from the un- 
derlying theory. Lattice simulations can, in principle, 
fix all of the coefficients. In order for us to investi- 
gate in some more detail the features of our potential 
and inspired by the PLM model mean-field type of ap- 
proximation we first assume 02 and bi to be positive 
and temperature independent constants while we model 
ai — a{T^ — T)/T, with T, a constant and a another 
positive constant. We will soon see that due to the inter- 
play between the hadronic states and i, T* need not to 
be the critical Yang-Mills temperature while oi displays 
the typical behavior of the mass square term related to 
a second order type of phase transition. 

The extrema are obtained by differentiating the poten- 
tial with respect to H and i: 
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dV 
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A. Small and Intermediate Temperatures 

At small temperatures the second term in Eq. ( p^ ) 
dominates and the only solution is £ = 0. A vanish- 
ing £ leads to a null P[£] yielding the expected minimum 
for H: 



e 



(15) 



Here £ and H decouple. 

We now study the solution near the critical tempera- 
ture for the deconfinement transition. For all the tem- 
peratures for which 



T^ai -I- Hbi = T^a{T^ -T) + Hbi > 0, 



(16) 



the solution for £ is still £ = yielding Eq. (|T^. The 
critical temperature is reached for 

T. = n + ^^. (17) 

The critical temperature can be determined via lattice 
simulations. Wc sec that within our framework the latter 
is related to the glueball (gluon-condensate) coupling to 
two Polyakov loops and it would be relevant to measure 
it on the lattice. At T = T^, £ = and 7J = A^/e. 
Let us now consider the case T = Tc + AT with 



AT 



< 1. 



(18) 



Expanding (£)^ at the leading order in AT/Tc yields: 
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(19) 



We used Eq. ( |l7| ) and Eq. (^ which relates the tem- 
perature dependence of H to the one of £. At high tem- 
peratures (see next subsection) (£) can be normalized to 
one by imposing a/2a2 = 1 and the previous expression 
reads: 



{£r = 



eaTj 



AT 

~t7 



4Tr - 3T, 



(1 - 2bi)Tc. 



AT 
-—.(20) 



For a given critical temperature consistency requires 61 
and T, to be such that: 



ATc - 3T* 



(1 - 2bi)T^ + 2biT, 



> 0. 



(21) 



The temperature dependence, in this regime, of the gluon 
condensate is: 



(i^) = — exp [~2b,{£)^] 



(22) 



Wc find the mean field exponent for £, i.e. £'^ increases 
linearly with the temperature near the phase transition 
p?! . Interestingly the gluon-condensates drops exponen- 
tially. The drop in the gluon-condensate is triggered by 
the rise of £ and it happens in our simple model exactly at 
the deconfining critical temperature. Although the drop 
might be sharp it is continuous in temperature and this 
is related to the fact that the phase transition is second 
order. Our findings strongly support the common picture 
according to which the drop of the gluon condensate sig- 
nals, in absence of quarks, deconfinement. 



B. High Temperature 

At very high temperatures the second term in Eq. (|lj) 
can be neglected and the minimum for t is: 




(23) 
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For H we have now: 
A* 



■ exp 



-26i 



2a2 



exp [-2bi] . (24) 



In the last step we normahzed [t] to unity at high tem- 
perature. In order for the previous solutions to be valid 
we need to operate in the following temperature regime: 



a 



(25) 



We find that at sufficiently high temperature (H) is ex- 
ponentially suppressed and the suppression rate is de- 
termined solely by the glueball - mixing term en- 
coded in V[£]. The coefiicient bi should be large (or 
increase with the temperature) since we expect a van- 
ishing gluon-condcnsate at asymptotically high tempera- 
tures. Clearly it is crucial to determine all of these coeffi- 
cients via first principle lattice simulations. The qualita- 
tive picture which emerges in our analysis is summarized 
in Fig. g. 
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FIG. 1: The thin line is the gluon condensate {H) normal- 
ized to A*/e as function of the temperature. The thick line 
represents the normalized to unity {£) as function of the tem- 
perature. We have chosen for illustration a = 1, bi = 1.45 
and Tc ~ 1.16A. 



IV. THE THREE COLOR THEORY 

Z3 is the global symmetry group for the three color 
case while ^ is a complex field. The functions V[£] and 
V\£] are: 



V[£] = a^\£\^ + a2\£\^ 
V[£] = bx\£\'' + 0{£^), 



a,{£' 



r^) + 



(26) 



with ai, 02, 03 and 61 unknown temperature dependent 
coefficients which can be determined using lattice data. 
In this paper we want to investigate the general relation 
between glueballs and £ so we will not try to find the best 
parameterization to fit the lattice data. In the spirit of 



the mean field theory we take 02, 03 and bi to be positive 
constants while ai = a{T^ - T) /T. With £ = |4e"^ the 
extrema are now obtained by differentiating the potential 
with respect to H , \£\ and (p: 



dV 
dH 
dV 
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^61 - SasKI cos(3</p) + 2a2|^r ) = 0, 



6|£|^ sin(3(p) = 0. 



(27) 



At small temperature the H /T^ term in the second equa- 
tion dominates and the solution is (|^|) = 0, {H) = A^/e 
and the last equation is verified for any {ip) , so we choose 
((ys) = 0. The second equation can have two more solu- 
tions: 



3a3_^ /_9_a3 
16 a| 



air - T*) biH 



4 02 



2ra2 



2a2T4 ' 



(28) 



whenever the square root is well defined (i.e. at suf- 
ficiently high temperatures). The negative sign corre- 
sponds to a relative maximum while the positive one to 
a relative minimum. We have then to evaluate the free 
energy value (i.e. the effective thermal potential) at the 
relative minimum and compare it with the one at € = 0. 
The temperature value for which the two minima have 
the same free energy is defined as the critical tempera- 
ture and is: 



41 
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(29) 



When 03 vanishes we recover the second order type criti- 
cal temperature Tc- To derive the previous expression we 
held fix the value of H to A^/e at the transition point. 
In a more refined treatment one should not make such 
an assumption. Below this temperature the minimum is 
still for {£) = and {H) = A^/e. 

Just above the critical temperature the fields jump to 
the new values: 



(Kl)--, {H) = -e^p[-2bM']- 



A^ 



Close but above Tc (i.e. T 
with 



03 
a2 



Tc + AT) we have: 
AT 



(30) 



(31) 



P 



aa2 



AkTc — 3T* 



as a2Tc - Abia{KTc - T,) ' 
aa2 + ai 



aa2 



(32) 



a positive function of the coefficients of the effective po- 
tential. In this regime 



A^ 

{H) = — exp 
e 



-2foi(— +p— ) 

a2 -Lc 



(33) 



5 



At high temperature we expect a behavior similar to the 
one presented for the two color theory case. A cartoon 
representing the behavior of the Polyakov loop and the 
gluon condensate is presented in Fig. ||. 

Since we are in the presence of a first order phase tran- 
sition higher order terms in Eq. ( p6|) may be important. 
However lattice simulations have shown that the behav- 
ior of the Polyakov loop, for 3 colors, resemble a weak 
first order transition (i.e. small a^) partially justifying 
our simple approach. 

The approximation for our coefficients is too crude and 
it would certainly be relevant to fit them to lattice sim- 
ulations. 

What we learnt is that the gluon condensate, although 
not a real order parameter, encodes the information of 
the underlying symmetry. More generally we have 
shown that once the map between hadronic states and 
the true order parameter is known we can use directly 
hadronic states to determine when the phase transition 
takes place and the order of the phase transition. For 
example by comparing Fig. |^ and Fig. ^ we immediately 
notice the distinct difference in the gluon condensate tem- 
perature dependence near the phase transition. 
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FIG. 2: A cartoon sketcliing the gluon condensate (H) nor- 
malized to A'^/e and the {\£\) (thick line) behaviors as function 
of the temperature. We have chosen for illustration as — 0.3, 
aa = 1, a = 2, 6i = 0.7 and Tc = 1.2A. 



V. CONCLUSIONS AND SELF CRITICISM 

Our simple model is able to account for many features 
inherent to the Yang-Mills deconfining phase transition. 
We related two very distinct and relevant sectors of the 
theory: the hadronic sector (the glueballs), and some 
dimensionless fields {£) charged under the discrete group 
Zn understood as the center of the underlying SU{N) 
Yang-Mills theory fill. 

The gluon-condensate is, strictly speaking, not an or- 
der parameter for the deconfining Yang-Mills phase tran- 
sition. However we have shown that the information en- 
coded in the true order parameter i is efficiently commu- 
nicated to the gluon condensate. Since the exponential 



drop of the condensate just above the Yang-Mills critical 
temperature is a direct consequence of the behavior of 
the true order parameter at the transition we can con- 
sider this drop as a strong signal of deconfincmcnt. This 
drop has already been used in various models for the 
Yang-Mills phase transition. We have also seen that the 
reduction in the gluon-condensate is associated to the 
increase of the Polyakov loop condensate i. The infor- 
mation about the order of the phase transition is also 
transferred to the behavior of the gluon condensate. In- 
deed from Fig. and Fig. ^ we deduce that the drop is 
continuous for the gluon condensate in the two color case 
while is discontinuous for the three color theory. Phys- 
ically the glueballs start decaying into gluons and this 
information is encoded in the term of the present 

theory. 

We now better understand the mechanism for trans- 
ferring information from the Yang-Mills order parameter 
to the physical states. 

It is important to stress that our model is very limited 
since we neglect the temperature corrections associated 
with the glueball gas as well as other dimensionless op- 
erators with different charges under Z^. Finally we did 
not include any of the excited glueball and pseudo-scalar 
glueball states present in the theory. Besides the temper- 
ature dependence of the coefficients in Vll] and V[i] is not 
known and we have just adopted the simplest guess con- 
sistent with mean-field theory. We also know that mean- 
field cannot be the whole story and corrections need to 
be included. 

It is worth mentioning that the Polyakov loop need not 
to be the only acceptable order parameter. For example 
using an abelian projection one can define a new (non 
local in the cromomagnetic variables) order parameter 
. Our model can be, in principle, modified to be able 
to couple the hadronic state to any reasonable Yang-Mills 
order parameter. 

The same holds true when considering the introduc- 
tion of quarks. Once identified a true order parameter 
for QCD with quarks we can first construct a model La- 
grangian which satisfies the ordinary symmetries at zero 
temperature for the hadronic states and then extended 
it to describe at the same time the order parameter and 
the hadronic states. Although i is not a good order pa- 
rameter when quarks are added to the theory since the 
Zn symmetry is explicitly broken we can still construct 
a theory containing i and the hadronic states (mesons 
and glueballs) provided we introduce explicit Zn break- 
ing terms. This approach might be relevant for under- 
standing RHIC physics Q |l9). 

Although the model is at a very early stage of devel- 
opment at the tree level some of the results are already 
fairly robust. For example the exponential drop of {H) 
as function of {£) is a prediction not expected to be very 
sensitive to different sources of corrections. We also note 
that the first order behavior of the deconfining three color 
Yang-Mills theory is directly inherited by the (H). In- 
deed this quantity is discontinuous at the phase transi- 
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tion for three colors while displays a smooth behavior in 
the two color case. We expect these results to be quite 
general. We also stress that they are connected to the 
saturation of trace anomaly and the Zj^ symmetry at 
the effective Lagrangian level when considering simulta- 
neously £ and H. 

By computing the temperature dependence of the coef- 
ficients in the present effective theory lattice simulations 
can test the validity of the present model while improving 
the present results. 
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